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lO I ground coupled to a Maxwell field via marginally stable modes. In the probe 

On I limit, we study the magnetic field effect on the holographic insulator/superconductor 

phase transition numerically and analytically. The condensate will be localized 
in a finite circular region for any finite constant magnetic field. Near the critical 
point, we find that there exists a simple relation among the critical chemical 
^ ' potential, magnetic field, the charge and mass of the scalar field. This rela- 

' tion indicates that the presence of the magnetic field causes the transition from 

insulator to superconductor to be difficult. 



We investigate the scalar perturbations in an anti-de Sitter soliton back- 



1 Introduction 

The AdS/CFT correspondence [1] provides a powerful theoretical method to understand 
the strongly coupled field theories in condensed matter physics. A holographic super- 
conductor (superfluid) model has been constructed recently in the work [2, 3] (for reviews 
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see [4]). The dual gravitational configurations are some anti-de Sitter (AdS) black holes 
with/without some charged matter contents. 

The phase transition studied in Refs.[2, 3] is actually a holographic superconduc- 
tor/metal phase transition. The holographic superconductor model can be simply con- 
structed by an Einstein-Maxwell theory with a negative cosmological constant coupled to 
a complex scalar field. In particular, when the temperature of the black hole is below a 
critical temperature, there are at least two distinct mechanisms leading to the black hole 
solution unstable to develop a scalar hair near the horizon [5]. And the condensation of 
the scalar hair induces the local U(l) symmetry breaking of the system, which gives the 
nonvanishing vacuum expectation value to the dual charged operator in the boundary field 
theory. Therefor, the U(l) symmetry breaking in gravity leads to a breaking of the global 
U(l) symmetry in the dual field theory. This results in a superconductor (superfiuid) phase 
transition. 

The holographic insulator/superconductor phase transition was first researched in Ref. [6]. 
A five-dimensional AdS soliton background [7] coupled to a Maxwell field and a scalar field 
was used to model the holographic insulator/superconductor phase transition at zero tem- 
perature. The normal phase in the AdS soliton is dual to a confined gauge theory with a 
mass gap which is reminiscent of an insulator phase [8]. As the chemical potential grows 
sufficiently up to a critical value, the instability is triggered, resulting in the emergence 
of the scalar hair which is dual to a superconducting phase in the boundary field the- 
ory. The holographic insulator /superconductor phase transition was also investigated in 
Refs.[9, 10, 11, 12, 13, 14, 15]. 

There are some discussions of holographic superconductors in the presence of magnetic 
field [5, 16, 17, 18, 19, 20, 21], but they mainly focused on superconductor/metal phase 
transition. Motivated by these studies and especially by our previous works [13, 22], we 
expect to explore how the magnetic field impacts the behavior of the condensate if a 
magnetic field is added into the AdS soliton spacetime. We are working in the probe limit, 
which means that the scalar field and Maxwell field have no back-reaction to the gravity 
background. 

In this paper, we first adopt the idea of marginally stable modes [2] to study the 
holographic insulator/superconductor phase transition in AdS soliton background at zero 
temperature. We employ the method introduced in Ref. [23] to study the quasinormal 
modes(QNMs) of the scalar perturbations of the system just like the one in Ref. [6]. At 
some critical values of the chemical potential fi and magnetic field B, the marginally 
stable modes will emerge. This represents that at the critical parameters, the AdS soliton 
background becomes unstable and will be fond of an AdS soliton background coupled with 
the hair of charged scalar fields. In particular, for some given mass and charge q, the 
special combination of the critical chemical potential fi and the critical magnetic field B 
satisfies a simple relation, i.e. q'^fi'^ — \ qB\ = A^, where A is a positive constant. The minus 
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sign in the combination is interesting because, as the magnetic field grows stronger, the 
critical chemical potential becomes higher. This indicates that the presence of the magnetic 
field causes the transition from insulator to superconductor to be difficult. Setting the 
magnetic field vanishing, the critical chemical potentials we derived are consistent with 
the results obtained by previous works [13, 22]. Actually, corresponding to various critical 
parameters, there are multiple marginally stable modes related to the nodes n = 1, 2, 3 ■ ■ ■ , 
which are unstable due to the oscillations of the field in the radial direction [2, 24]. Taking 
advantage of the shooting method, we plot the profile of the scalar field depending on the 
radial direction. From these diagrams, one can intuitively see the "nodes" of the scalar field. 
Furthermore, using the variational method for the Sturm-Liouville eigenvalue problem [25], 
we analytically study the holographic insulator/ superconductors phase transition, following 
the previous work in Ref. [13]. Near the critical point, we show a simple relation among 
the critical chemical potential, magnetic field, the charge and mass of the scalar field. 

The paper is organized as follows: in Sec. 2, we introduce the AdS soliton background 
and obtain the equations of motion in the probe limit. We study the marginally stable 
modes in Sec. 3. In Sec. 4 the system is solved by shooting method. In Sec. 5 we extract the 
relation of the four parameters at the critical phase transition points using the variational 
method for the Sturm-Liouville eigenvalue problem. Conclusions and discussions are drawn 
in Sec. 6. 



2 The Background 

We construct the model of the holographic insulator/superconductor phase transition with 
the Einstein-Maxwell-scalar action in five-dimensional spacetime: 

S = J d'xV^iR + ^ - - |V^^ - tqA,ij\' - m^j^H, (1) 

where L is the radius of AdS spacetime. When the Maxwell field and scalar field are absent, 
the above action admits the AdS soliton solution [7] : 

ds^ = L^j^ + r\-dt^ + dp" + p'de^) + f{r)dx^. (2) 

where /(r) = — rg/r^.^ The asymptotical geometry approaches to i?^'^ x near the 
boundary. Moreover, the Scherk-Schwarz compactification x ~ X + '^L/tq is needed to 
obtain a smooth geometry. This gives a dual picture of a three-dimensional field theory 
with a mass gap, which resembles an insulator in the condensed matter physics. The 
geometry in (r, x) directions just looks like a cigar whose tip is at r = tq. The temperature 
in this background is zero. 

^ We will work in polar coordinates dx^ + dy^ = dp^ + p^dO^ in this paper. 



3 



The equations of motion (EoMs) of matter fields are 



( - iqA^) ( - iqA'')^ - rr?^ = 0, (3) 
V^F''^ = iq[ip{V - iqA^')ip - ^{\/^ + iqA'')ip]. (4) 

The boundary conditions for the matter fields near the infinity r — )■ oo of the AdS soliton 
are 



A = /.-^ + ..., (6) 

where = i = 1,2, and 0(j) are the corresponding dual operators of in the 

boundary field theory. The conformal dimensions of the operators are A± = 

fj, and Pe are the corresponding chemical potential and charge density in the boundary field 

theory. 

Following Ref. [6], we work in the probe limit, where the Maxwell field and scalar field 
do not back react on the background metric. When there is no condensate, i.e. ip = 0, 
EoMs have the solution 

A = ndt+ ^Bp^de. (7) 

In addition to a constant chemical potential p, adding a constant magnetic field B to the 
Maxwell field is also consistent in the probe limit. 

As the scalar field ip does not vanish, we must solve the coupled EoMs. However, 
near the critical point of phase transition the scalar field is nearly zero, we can treat ip 
as a probe into this background which is a neutral AdS soliton with a constant electric 
potential and magnetic field. We are interested in axisymmetric solutions in which all 
fields are independent of 6. We consider an Ansatz of the form ip = F{t,r)H{x)U{p). 
Substituting ip into Eq.(3) and making the separation of the variables, we can obtain 

d^F{t,r) ^ ^3 ^drf^dF{t,r) L'^ d'^F{t,r) ^ 2iqpL'^ dF{t,r) 



Qy,2 r f dr fr"^ dt'^ fr"^ dt 

- mh' - ^ - k')F{t, r) = 0, (8) 

where and k"^ are the eigenvalues of the following equations, respectively: 



yHix). (9) 



^^[4^(^)1 -i^^^V^(p)^-^^^(p), (10) 
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where A = 2rol/L, I G Z owing to the periodicity of H{x) = H{x + t^L/tq). U{p) solves 
the equation for a two-dimensional harmonic oscillator with frequency determined by B 
and k"^ = n\qB\, n G Z"*". We expect that the lowest mode / = 0,n = 1 will be the first to 
condense and result in the most stable solution after condensing. We can also set L = 1 
and ro = 1 without loss of generality. We finally obtain the equation of motion of F{t,r): 

d^F{t,r) ^ ^3 ^ drf ^ dF{t,r) _ d^F{t,r) ^ 2iqfiL^ dF{t, r) 
dr'^ r f dr fr"^ df^ fr"^ dt 

+Tt(5V -\qB\- m^r^)F{t, r) = 0. (11) 

Note that in this case, U{p) = exp( "^^^^^ ), so for any finite magnetic field, the super- 
conducting condensate will be localized to a finite circular region. As the magnetic field 
becomes smaller, the region grows until it occupies the whole plane, which can be seen in 
the profile of U{p) by setting i? — > 0. We can also see similar phenomena in other gravity 
background [5, 17]. 



3 Critical behavior via Quasinormal Modes 

To reveal the stability of a spacetime background, an effective method is to analyze the 
QNMs of the perturbations in the fixed background (for reviews, see Refs.[26, 27, 28]. The 
temporal part of the QNMs behave like e"*"^*. Therefore, if the imaginary part of the QNMs 
is negative, the mode will decay in time and the perturbation will ultimately fade away, 
indicating that the background is stable against this perturbation. In sharp contrast, if 
the imaginary part is positive, the background is unstable. The critical case is that if the 
perturbation has a marginally stable mode, i.e. u = 0, one always expects that this is a 
signal of instability, where a phase transition may occur [2]. 

In order to study the phase transitions in the background, we further define F{t,r) = 
e~*'^*i?(r), and Eq.(ll) becomes 

R "(r) + (^ + '(r) + ^[ico + qpf - \qB\ - mV]i?(r) = 0, (12) 

where a prime denotes the derivative with respect to r. 

Taking advantage of the Horowitz and Hubeny's method [23] to study these QNMs, we 
find it is convenient to work in the ^-coordinate where z = \/r. In this new coordinate, 
the infinite boundary is now at 2; = 0, while the tip is ai z = Zq = 1 / Tq = 1 . Equation (12) 
becomes 

R "{^) + [-- + 4rT]^ '(^) + -Ar-yi^ + ^^'f - ^'\<iB\ - ^']R{^) = 0, (13) 
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where we make the prime denote the derivative with respect to z from now on. Following 
the procedure of Horowitz and Hubeny [23], we multiply z^ji^z^ji^z — 1) to both sides of 
the above equation, and we obtain 

S{z)R '\z) + \z) + = 0, (14) 

2 — 1 [z — vy 

where the coefficients are given by 

S{z) = (15) 
z — 1 

T{z) = -z^f{z)+z^f'{z), (16) 



V{z) 



z^{u + q^f - z^\qB\ - nf 



[z-1). (17) 



S{z),T{z) and V{z) are all polynomials and can be expanded to a finite order as 

J 

S{z) = J2s,{z-iy, (18) 
i=o 

where J is a finite integer. T{z) and V{z) can be polynomially expanded similarly. 

Because of the absence of a black hole horizon in the AdS soliton, the boundary condi- 
tion at the tip is a finite quantity. This motivates us to find a solution like 



R(z) = lim y a„(z- 1)". (19) 

n=0 

Substituting Eqs.(19) and (18) into Eq.(14) and comparing the coefficients of {z — 1)" for 
the same order, we find that 

^ n— 1 

an = -—^[k{k-l)sn-k + ktn-k + Vn-k\ak, (20) 

Hn = n(n — l)so + nto + vq = —Arf . (21) 

We set ao = 1 for simplicity because of the linearity of Eq.(14). The boundary condition 
for the scalar field at 2; = is 



N 

i?(0)= hm J]a„,(-1)" = 0. (22) 



n=0 



And the algebraic equation (22) can solve the modes u. 

We restrict q = 1 just like the one in Ref. [6] in the numerical calculations. In practice, 
we will expand R{z) to a large order which is = 300 in this paper. Because of the 
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fact that the real part of the frequency indicates the energy of the mode, we only take 
care of the QNMs with a positive real part. From the numerical results which are plotted 
in Figure. (1), we find that when a marginally stable mode arises, the square of chemical 
potential /i^ and magnetic field B satisfy a simple linear relation, whose slope is particularly 
one, and the intercept just gives the square of critical chemical potential in the absence 
of magnetic field, which is consistent with our previous work [22]. For convenience, we 
define n"^ — \B\ = A"^ which determines the stability of the system. The minus sign in the 
combination is interesting, which means as the magnetic field grows stronger, the critical 
chemical potential becomes higher. This indicates that the presence of the magnetic field 
causes the transition from insulator to superconductor to be difficult. By increasing A to 
some critical value, the marginally stable modes will emerge, or rather, a phase transition 
may occur. 

Table 1: The first three lowest-lying A„'s for various mass squares obtained from the 
calculation of the marginally stable modes. 
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(B.) 



Figure 1: (A.) The square of the critical chemical potentials for the first (red dots) and the 
third (purple quadrate points) lowest-lying modes versus magnetic field at = —15/4. 
(B.) The first three lowest-lying A^'s for the marginally stable modes versus of the 
scalar field. The points are numerically obtained from the QNMs. 



Table (1) shows the first three lowest-lying critical A„'s in which the index n denotes 
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the "overtone number" for various m^'s. Marginally stable modes corresponding to higher 
overtone numbers n = 1, 2 ■ ■ ■ may also appear. However, they are unstable due to the 
oscillations in the z direction, which can be explicitly seen in the next section. The nodes 
n = 0, 1, 2 can also be intuitively seen in the next section. The part B of Fig.(l) shows A„'s 
of the marginally stable modes for various squared mass of the scalar field. Turning off 
the magnetic field, we recover the results previously derived in Ref.[22]. Therefore, despite 
that we do not exactly know the phase structures through the marginally stable modes, 
they can reveal the onset of the phase transition in practice. 



4 Critical behavior via the shooting method 

Shooting method is an alternative way to study the critical behavior of the phase transition, 
which has been widely used in the previous studies on holographic superconductors [4]. In 
this section, we will make use of the shooting method to study the critical behavior in AdS 
soliton background, especially to plot the profile of the scaler field. We will also compare 
it with the above quasinormal modes method. 

We focus on the static case in which F is independent of t; we can get equation of 
motion by simply setting a; = in Eq.(13): 

R "{z) + [-1 + '{z) + ;r^[^'(9V - \qB\) - m']R{z) = 0. (23) 

Near the boundary R{z) behaves as 

i?|,^0 = <^(l)z2-vW + ^(2)^2+vW ^ _ _ _ _ ^24) 

It is well known that when < ■\/4 + m? < 1, the scalar field admits two different 
quantizations related by a Legendre transform [29]. <y5*^*^ can either be interpreted as a 
source or an expectation value. In this paper, we will only consider the case where the 
faster falloff is dual to the expectation value, i.e. cp^^"^ = 0. 

We now study the behavior of the solution near the tip z = 1. Then Eq.(23) becomes 

R \z) - ^R \z) + ^(X^^(^) = 0, (25) 

where k = g^/i^ — \qB\ — rn? ^ 0.^ It has a solution near 2: = 1 as 

Rl^i = a + l3 log( v^(l - z)) + . . . , (26) 

where a and /3 are two constants. Since we want the field to be finite at the tip, we choose 
the boundary condition by setting /3 = 0. 

^Note that we neglect the particular case k = 0. One can also consider this case, but it does not matter. 
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To make use of the shooting method, we begin with an initial value of R{z) at the 
tip Zq and then calculate the EoM of R{z) [Eq.(23)] numerically provided that the infinite 
boundary condition Lf^^^ = is satisfied. We also take q = 1 for simplicity. We can set 
i?(l) as an arbitrary constant due to the linearity of Eq.(23). In addition, near the critical 
point of the phase transition the quantity of R{z) is very close to zero, therefore, we can 
impose the following initial conditions at the tip z = 1: 

i?(l) = 0.001, i? '(1) = (27) 

For a given m^, only for certain values of fi and B do we get to satisfy the boundary 
conditions. 

Note that the chemical potential fi and magnetic field B appear only as a whole g^/i^ — 
\qB\ in Eq.(23) and especially in n. And inspired by the result from marginally stable 
modes that only the combination A"^ = n"^ — \B\ determines the occurrence of instability, 
we expect that this result may be confirmed by our shooting method. Just as expected, we 
find that only the appropriate can trigger the condensate. Furthermore, the values of 
Ac obtained from the shooting method are perfectly consistent with the A.„'s in Table (1) 
derived from the "marginally stable modes" method. So we will not distinguish between 
them from now on. 

Figure. (2) exhibits the multiple marginally stable curves of the scalar fields for various 
m^. For example, the first three lowest-lying modes in the plot of = —15/4 are in the 
sequence Aq < Ai < A2. The red line is dual to the minimal value of A^, which has no 
intersecting points with the R{z) = axis at nonvanishing z. Therefore, we consider the 
mode corresponding to Ac ~ 1.8816 as a mode of node n = 0. Further, we regard the green 
line dual to Ac ~ 4.2341 and blue line related to Ac ~ 6.6156 as modes with nodes n = 1 
and n = 2, respectively. However, the green and blue lines are thought to be unstable 
due to the fact that the radial oscillations in 2;-direction of R{z) will cost energy [24]. In 
addition, the above discussions also hold for other diagrams in Fig. (2). It is interesting to 
note that the marginally stable curves corresponding to different Ac's are just the division 
of the curves of R{z) with different nodes. 

5 Critical behavior via the Sturm-Liouville method 

From the above numerical analysis, we can see that when the combination of chemical 
potential /i and magnetic field B, which is A^ = /i^ — \B\, exceeds a critical value Aq for 
given mass, the condensations of the operators will turn out. This can be regarded as a 
superconductor (superfluid) phase. However, when less than Aq, the scalar field is vanishing 
and this can be interpreted as the insulator phase because this system has a mass gap, 
which is due to the confinement in the (2+l)-dimensional gauge theory via the Scherk- 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

z z 

Figure 2: The marginally stable curves of the scalar field corresponding to various critical 
Ac's in the cases of different mass square. The critical Ac's for different curves are Aq < 
Ai < A2 in the sequence . 

Schwarz compactificaiton. Therefore, the critical parameters satisfied, Aq = /i^ — are 
the turning points of this holographic insulator/superconductor phase transition. 

In this section, using the variational method for the Sturm-Liouville eigenvalue problem 
[25], we analytically study the phase transition just following the procedure in Ref.[13]. 
Especially, we look forward to giving an explicit demonstration to our previous numerical 
results and trying to find an approximate function to relate parameters {q, fi, B,m^} at 
the critical phase transition point. We start from Eq.(23) and also focus on the case 
(f^^^ = 0. The operator 0(2) is normalizable when > = —4, where m'^p is the 
Breitenlohner-Freedman bound of the mass square of scalar field in the AdS spacetime. 

Following the steps in Ref.[13], we introduce a trial function T{z) into R{z) near z = 
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as 



^ (0(2))^^+^^+'"^r(^). (28) 

The boundary conditions for T{z) are r(0) = 1 and r'(0) = 0. It is easy to obtain the 
EoM of rf^) as 



2:^(8 + 771^ + 4 a/4 + m2) 



+ ^ l"^-^^-v-^-^ r(.) = 0. (29) 



The EoM of T{z) can be rewritten as 
d 



r 1+2^4+^(^4 _ ^) p /.X. ^ ^3+2VI+^^g + ^2 ^ 4^4+^) r(^) 

a2 ^ ^ ' ^ „, ' 



K 



-2^+2^'^(gV - \qB\)T{z) = 0. (30) 
Q 

The minimum eigenvalue of q'^fi'^ — \qB\ can be obtained by varying the following functional: 

2 2,, [Uz{Kr'^ + PT^) 

?V-k5| = ^° \ ^=e(7,m2). (31) 

In order to estimate it, we try to set T{z) = 1 — 72^. Thus, we obtain a particular function 
g^/i^ — \qB\ = ^(7,m^), from which we can get the minimal value of g^/i^ — \qB\ for given 
mass. It will be subtle if the function ^(7;''^^) has more than one minimum. Fortunately, 
we can see from part A of Fig. (3) that there indeed exists only one minimum for given 
mass square. We denote the minimum for given mass square as A[m^]. 



Table 2: The critical parameters g^/i^ — \qB\ for various mass squares obtained from the 
QNMs, shooting method, and S-L method. 





2 

m = 


-15/4 


= -3 


= -7/3 


= —1 


QNMs 


11 


^849 


2.3963 


2.6903 


3.1346 


Shooting 


l.t 


^816 


2.3964 


2.6904 


3.1347 


S-L 


l.J 


^904 


2.3986 


2.6928 


3.1373 



In Table (2), we list the critical parameters Aq for various mass squares obtained from 
the three methods, i.e. the calculation of the marginally stable modes, the shooting 



11 




Figure 3: (A.) ,^ as a function of 7 and m^. (B.) The critical chemical potenial |g/i| versus 
\qB\ and m?. (C.) The charge q versus mass square at fixed chemical potential /i = 1 
and magnetic field B = 0.5. (D.) The critical chemical potential /i versus the charge q at 
fixed mass square = — 1 and magnetic field B = 0.5. 

method, and the Sturm-Liouville (S-L) method. The analytical results are in good agree- 
ment with numerical values. So we can trust the analytical treatment, especially the choice 
of r{z) = 1 - 7^2. 

We finally obtain the simple relation for parameters q,^,B,m? at the critical phase 
transition point as 

gV_ = A[m2]. (32) 

This simple relation can explain the main results we obtained. We plot the function 
= \qix\{\qB\,im?) in part B of Fig.(3). The relation can be regarded as a function of 
three arbitrary independent variables in {g, i?,m^}. We extract the relation between 
the charge q and mass square m? for the systems whose phase transitions take place at 
fixed chemical potential and magnetic field B. Part C of Fig. (3) shows schematically 
that when the mass square of the scalar field grows, the charge must grow too. Similarly, 
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it can be seen from part D of Fig. (3) that the critical chemical potentials are sufficiently 
depressed by the charge of the scalar field for fixed mass and magnetic field. 

6 Discussions and Conclusions 

In the probe limit, we studied the holographic insulator /superconductor transition in an 
external constant magnetic field both numerically and analytically. When the magnetic 
field is absent, the condensate fills the plane homogeneously. On the other hand, noting 
that the profile of scalar field 'ijj{z, p) = R{z) exp( ~^'^^^^ ) where the profile of R{z) in 
superconducting phase is shown in Fig. (2), we see that for any nonvanishing magnetic 
field, the superconducting condensate will be localized to a finite circular region. As the 
magnetic field becomes smaller the region grows until it occupies the whole plane, which 
can be seen in profile of U{p) by setting 5—7-0. 

Despite that we do not exactly know the phase structure by investigating the marginally 
stable modes of the scalar field, they can actually indicate the onset of the phase transi- 
tion, which means that the neutral AdS soliton will become unstable to develop charged 
scalar hairs in this AdS soliton background when the parameters are beyond the critical 
values. Our results show that only the parameter = g^/x^ — \qB\ triggers the instability. 
Making use of the shooting method to numerically plot the behavior of the scalar field in 
the radial direction, one can intuitively see the nodes of the marginally stable modes. Fur- 
thermore, taking advantage of analytical approach we directly obtained a simple relation 
which constrains parameters {g,/i, -B,m^} at the critical point, i.e. q^p? — \qB\ = A[m^]. 
This relation is the main result of our paper, and it can be used to explain many properties 
of the model. 

Note that in this paper, we have limited ourselves to the choice of ^/''•^^ = and A = 0, 
and n = 1. It is interesting to study the cases with vanishing -ip^^^ and other (A,n)'s. Of 
course, it is also quite significant to investigate the back reaction of matter fields and to 
draw a full phase diagram of the holographic insultor /superconductor transition. These 
discussions in the paper are desired to extend to the background of modified gravity, such 
as including Chern-Simons term and term[30, 31]. We expect to report on further 
studies on these and relevant issues. 
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